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A rigorous three-dimensional relativistic equation for quark-antiquark bound states at finite 
temperature is derived from the thermal QCD generating functional which is formulated in the 
coherent-state representation. The generating functional is derived newly and given a correct path- 
integral expression. The perturbative expansion of the generating functional is specifically given by 
means of the stationary-phase method. Especially, the interaction kernel in the three-dimensional 
equation is derived by virtue of the equations of motion satisfied by some quark-antiquark Green 
functions and given a closed form which is expressed in terms of only a few types of Green functions. 
This kernel is much suitable to use for exploring the deconfinement of quarks. To demonstrate the 
applicability of the equation derived , the one-gluon exchange kernel is derived and described in 
detail. 
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I. INTRODUCTION 

Quantum Chromodynamics (QCD), as a strong interaction theory of quarks and gluons, has a distinctive property 
' of asymptotic freedom and infrared slavery which makes the quarks and gluons to be confined in hadrons. It is 
, widely believed that in extreme conditions, i.e., at high temperature and/or high density, the quarks and gluons 
\^ ' would be deconfined from hadrons and form a new matter, the quark-gluon plasma (QGP). It is highly expected 
■f—i ' that the QCD phase transition from hadrons to QGP would take place and be observed in the high energy heavy ion 
' collisions at RHIC [1-3]. Theoretically, to predict the QCD phase transition, many efforts have been made by using 
[ different approaches such as the lattice simulation, the effective field theory, the hydrodynamic model and etc. [4-14]. 
1^ . According to the prediction of the lattice QCD calculations, the phase transition could occur when the colliding system 
' reaches the temperatures 150-170 MeV [2]. Since the quarks and gluons are confined in hadrons which exist as the 
. bound states of quarks and/or gluons, it is obvious that a proper approach of investigating the quark deconfinement 
"*7^ ' is to start from an exact relativistic equation for quark and/or antiquark bound states at finite temperature. In this 
Oh, paper, we are devoted to deriving a rigorous three-dimensional relativistic equation of Dirac-Schrodinger type for 
^ quark-antiquark (qq) bound states at finite temperature. In particular, the interaction kernel in the equation will be 
given a closed and explicit expression which will be derived by following the procedure described in Refs. [15-17]. 
This constitutes the main purpose of this paper. Clearly, if the interaction kernel and the equation can be calculated 
by a suitable nonperturbation method, one can exactly determine at which temperature the quark and antiquark will 
^ , be deconfined from mesons. 
d In this paper, we intend to derive the aforementioned equation and kernel from the thermal QCD generating 

functional formulated in the coherent-state representation. For this derivation, we need first to give a correct expression 
of the generating functional in the coherent-state representation. This constitutes another purpose of this paper. Here 
it is necessary to mention that the corresponding path integral expressions for the partition functions and generating 
functionals given in the previous literature are not correct [18-22]. The incorrectness is due to that in the previous path 
integral expressions, the integral representing the trace is not separated out on the one hand and the time-dependence 
of the integrand in the remaining part of the path integral is given incorrectly on the other hand. Such path integral 
expressions can only be viewed as a formal symbolism because in practical calculations, one has to return to the 
original discretized forms which lead to the path integral expressions. If one tries to perform an analytical calculation 
of the path integrals by employing the general methods and formulas of computing functional integrals, one would 
get a wrong result. The partition functions and generating functionals for many-body systems discussed in quantum 
statistics were rederived in the coherent-state representation and given correct functional-integral expressions in the 
author's previous paper [23]. These expressions are consistent with the corresponding coherent-state representations of 
the transition amplitude and the generating functional in the zero-temperature quantum theory [24-26] . Particularly, 
when the functional integrals are of Gaussian type, the partition functions and the generating functionals can exactly 
be calculated by means of the stationary-phase method [25-27]. For the case of interacting systems, the partition 
functions and finite-temperature Green functions can be conveniently calculated from the generating functionals by 
the perturbation method. The coherent-state representation of the partition function and the generating functional 
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given in quantum statistics is now extended to the thermal QCD in this paper, giving a correct formulation for the 
quantization of the thermal QCD in the coherent-state representation. 

The remainder of this paper is arranged as follows. In Sect. 2, we quote the main results given in our previous paper 
for the quantum statistical mechanics. These results may straightforwardly be extended to the quantum field theory. 
In Sect. 3, we describe the coherent-state representation of the thermal QCD in the first order (or say, Hamiltonian) 
formalism. In Sect. 4, the quantization of the thermal QCD is performed in the coherent-state representation 
by writing out explicdtly the generating functional of thermal Green functions. To demonstrate the applicability and 
correctness of the generating functional, we pay our attention to deriving the pcrturbative expansion of the generating 
functional in the coherent-state representation. Section 5 will be used to establish the three-dimensional equation 
obeyed by the qq bound states at finite temperature. Section 6 serves to derive the closed expression of the interaction 
kernel appearing in the three-dimensional equation. In Sect. 7, the one-gluon exchange kernel and Hamiltonian will 
be discussed in detail. In the last section, some concluding remarks will be made. In Appendix, the pcrturbative 
expansion of the generating functional given in Sect. 4 will be transformed to the corresponding one represented in 
the position space. 



II. PATH INTEGRAL FORMULATION OF QUANTUM STATISTICS IN THE COHERENT-STATE 

REPRESENTATION 

First, we start from the partition function for a grand canonical ensemble which usually is written in the form 
[19-22] 

Z = Tre-f^^ (1) 
where (3 = -j^ with k and T being the Boltzmann constant and the temperature and 

K = H -liN (2) 

here /x is the chemical potential, H and A'' are the Hamiltonian and particle-number operators respectively. In the 
coherent-state representation, the trace in Eq. (1) will be represented by an integral over the coherent states. To 
determine the concrete form of the integral, for simplicity, let us start from an one- dimensional system. Its partition 
function given in the particle-number representation is 

oo 

Z=J2{n\e-^^\n). (3) 

Then, we use the completeness relation of the coherent states [18-28] 

j D{a*a) \a><a* \=1 (4) 

where | a > denotes a normalized coherent state, i.e., the eigenstate of the annihilation operator a with a complex 
eigenvalue a [18-28] 

a\a) = a \a) (5) 

whose Hermitian conjugate is 

{a*\a^ = a* {a*\ (6) 

and D{a*a) symbolizes the integration measure defined by [18-28] 

j-^, » ) — f j^da*da, for bosons; . ^ 

^ ' ^ da* da, for fermions. ^ ' 

In the above, we have used the eigenvalues a and a* to designate the eigenstates \a) and (a*|, respectively. It 
is emphasized that since we use the normalized cigcnfunction of the coherent state whose expression in its own 
representation will be shown in Eq. (15), the completeness relation in Eq. (4) has the ordinary form as we are 
familiar with in quantum mechanics. Inserting Eq. (4) into Eq. (3), we have 
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Z = J2 D{a*a)D{a'*a') {n \ a') {a'*\e-'^^ \a) {a* \ n) (8) 

where 

(n I a') = -^{aTe-'^"'^' ^ > 

are the energy eigenfunctions given in the coherent-state representation (Note: for fermions, n = 0, 1) [20-26]. The 

both eigenfunctions commute with the matrix element ^o'* | e~^-^ \ because the operator K(a'^,a) generally is 

a polynomial of the operator o+a for fermion systems. In view of the expressions in Eq. (9) and the commutation 
relation [20- 26] 

a*a' = ±a'a* (10) 

where the signs " +" and " — " are attributed to bosons and fermions respectively, it is easy to see 

(n I a') {a* \ n) = {±a* \ n){n \ a') . (11) 

Substituting Eq. (11) in Eq. (8) and applying the completeness relations for the particle-number states and coherent 
ones, one may find 

Z = J D{a*a){±a*\e-'^^\a) (12) 

where the plus and minus signs in front of a* belong to bosons and fermions respectively. 

To evaluate the matrix element in Eq. (12), we may, as usual, divide the "time" interval [0, /3] into n equal and 
infinitesimal parts, /? = ne. and then insert a completeness relation shown in Eq. (4) at each dividing point. In this 

way, Eq. (12) may be represented as [19-27] 



n-l 

Z = jD{a*a) n D{a*ai) {±a*\e-'^ \an-i) {al_,\e-'^ \an-2) ■ 

i=l 



(13) 



X «+i| e-^^ \ai) (all e-^^ |ai_i) • • • (a^j e"^^ \a) 
Since e is infinitesimal, we can write 

g-e^(a+,a) _ ^ _ sK(a+,a) (14) 

where K{a'^,a) is assumed to be normally ordered. Noticing this fact, when applying the equations (5) and (6) and 
the inner product of two coherent states [19-27] 

(a* I tti-i) = e-5<«'-i<-i'»'-i+<»^-i (15) 

which suits to the both of bosons and fermions, one can get from Eq. (13) that 

n— 1 n 

Z = /D(a*a)e-''*"/ J] D{a*a,) exp{-e K{alai-i) 

(16) 

i=l i=l 

where we have set 

±a* = a* , a — Go- (17) 

It is noted that the factor e~" " in the first integrand comes from the matrix elements (±a*| a„_i) and {al \ a) and the 
last sum in the above exponent is obtained by summing up the common terms —^a*ai and — |a*_jaj_i appearing in 
the exponents of the matrix clement (a* | aj_i) and its adjacent ones {a*^^ \ ai) and {a*_i \ ai-2)- As will be seen in 
Eq. (21), such a summation is essential to give a correct time-dependence of the functional integrand in the partition 
function. The last two sums in Eq. (16) can be rewritten in the form 
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n n—1 

i=i i=i (-^g^ 



= 2a„a„_i + ^a^ao + ^ 1. J^i - ai ( — 



Upon substituting Eq. (18) in Eq. (16) and taking the limit e — > 0, we obtain the path-integral expression of the 
partition functions as follows: 

Z = y£)(a*a)e-"*"yj)(a*a)e^(«*'") (19) 

where 

Yl ^da* {T)da{T), for bosons; 

J)(a*a) = { Vr''^ */ / N r . • (20) 
^ ^ Y[da {T)da[T), for fermions ^ ' 

T 

and 

I{a*,a) = ia*(/3)a(/3) + ia*(0)a(0) - dT[la*{T)a{T) 

^^a*{T)a{r)+K{a*{T),a{r))] (21) 
= a* (/3)a(/3) - dr[a* {T)a{T) + K{a* (r), a(T)) 

where the last equality is obtained from the first one by a partial integration. In accordance with the definition given 
in Eq. (17), we see, the path-integral is subject to the following boundary conditions 

a{(3) = ±a,a{Q) = a (22) 

where the signs "+" and "— " are written respectively for bosons and fermions. Here it is noted that Eq. (22) does 
not implies a(/3) = ±a and a*(0) = a* . Actually, we have no such boundary conditions. 

For the systems with many degrees of freedom, the functional-integral representation of the partition functions may 
directly be written out from the results given in Eqs. (19) -(22) as long as the eigenvalues a and a* are understood 
as column matrices a = {ai,a2, ■ ■ ■ ,ak,- ■ ■) and a* = (a*, ■ • •, a^, ■ • •). Written explicitly, we have 



(23) 



where 



Z = J D{a*a)e-<^' J ^{a* a)e^^''' '"^ 

n ^da^dak , for bosons; 
^) { Ylda^dak , for fermions, '■^^^ 

k 

Yl ^dal{T)dak{T) , for bosons; 
'^{a a) — { Yldal{T)dak{T) , for fermions '-^^^ 

fer 

I{a*, a) = 4(/3)afc(/3) - / dT[aUT)dk{T) + Jr(4(r), afc(r))]. (26) 
Jo 

The boundary conditions in Eq. (22) now become 

aUl3) = ±al , ak{0) = ak. (27) 

In Eqs. (23) and (26), the repeated indices imply the summations over k. If the k stands for a continuous index as 
in the case of quantum field theory, the summations will be replaced by integrations over k. 

It should be pointed out that in the previous derivation of the cohcrcnt-statc representation of the partition func- 
tions, the authors did not use the expressions given in Eqs. (16) and (18). Instead, the matrix element in Eq. (15) 
was directly chosen to be the starting point and recast in the form [18-22] 



and 
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K I a,.,) = exp{-|[a|(^^^) - (^^-^)a,_i]}. (28) 
Substituting the above expression into Eq. (13) and taking the hmit e ^ 0, it fohows [18-22] 

Z = J ^{a*a) exp{- dri^a* {tMt) - la*{T)a{T) + K{a%T), a(T))]}. (29) 

Clearly, in the above derivation, the common terms appearing in the exponents of adjacent matrix elements were not 
combined together. As a result, the time-dependence of the integrand in Eq. (29) could not be given correctly. In 
comparison with the previous result shown in Eq. (29), the expression written in Eqs. (19)-(21) has two fimctional 
integrals. The first integral which represents the trace in Eq. (1) is absent in Eq. (29). The second integral is defined 
as the same as the integral in Eq. (29); but the integrand are different from each other. In Eq. (19), there occur two 
additional factors in the integrand : one is ° which comes from the initial and final states in Eq. (13), another 
is e2[''*(^)"('^)+"*(°)"('')] in which a*(/3) and a(0) are related to the boundary conditions shown in Eq. (22). These 
additional factors are also absent in Eq. (29). As will be seen soon later, the occurrence of these factors in the 
functional-integral expression is essential to give correct calculated results. 

To demonstrate the correctness of the expression given in Eqs. (23)-(27), let us compute the partition function for 
the system whose Hamiltonian is of harmonic oscillator-type as we meet in the cases of ideal gases and free fields. In 
this case, 

K{a*a) = w/jaj^afe (30) 
where Wfe = £fc — /U with Sk being the particle energy and therefore Eq. (26) becomes 

I{a*,a) = al{p)au{(3) - / dT[al{T)ak{T) + W;t4(r)afe(r)]. (31) 
Jo 

By the stationary-phase method which is established based on the property of the Gaussian integral that the integral 
is equal to the extremum of the integrand which is an exponential function [25-27] , we may write 

I)(a*a)e^("*'") = e^°(°*'") (32) 

where /o(a*, a) is obtained from /(a*, a) by replacing the variables a*j^{T) and afe(r) in J(a*, a) with those values which 
are determined from the stationary condition 61 {a*, a) = 0. From this condition and the boundary conditions in Eq. 
(27) which implies ^a^(/3) = and 6ak{0) = 0, it is easy to derive the following equations of motion [24-26] 

dfc(T) + Wfeafe(T) = 0, afe(r) - WfeaJ(T) = 0. (33) 

Their solutions satisfying the boundary condition are 

a,(r) = ake-'^^- , 4(t) = ia^e'^-^^"'^). (34) 

On substituting the above solutions into Eq. (31), we obtain 

Io{a*,a) = ±4afee-'^'=^ (35) 

With the functional integral given in Eqs. (32) and (35), the partition functions in Eq. (23) become 

7 - { I £'(a*a)e-"^""'(i-'="''"") , for bosons; , . 

° ~ V i?(a*a)e-«^«^(i+«"''"'=) , for fermions. ^ ' 

For the boson case, the above integral can directly be calculated by employing the integration formula [18]: 

J i?(a*a)e-«*(^"-*)/(a) = jf{X-'b) (37) 
The result is well-known, as shown in the following [20-22,29] 
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^0 = n ^''^ 

For the fermion case, by using the property of Grassmann algebra and the integration formulas [21-26]: 

j da = j da* = , j da*a* = j daa = 1 (39) 

it is easy to compute the integral in Eq. (36) and get the familiar result [21-23,29] 

Zo = l[{l + e-^"'') (40) 

k 

It is noted that if the stationary-phase method is applied to the functional integral in Eq. (29), one could not get the 
results as written in Eqs. (38) and (40), showing the incorrectness of the previous functional-integral representation 
for the partition functions. 

Now let us turn to discuss the general case where the Hamiltonian can be split into a free part and an interaction 
part. Correspondingly, we can write 

K{a*,a) = Ko{a*,a)+Hi{a*,a) (41) 

where Ko{a* ,a) is the same as given in Eq. (30) and Hi{a*, a) is the interaction Hamiltonian. In this case, to evaluate 
the partition function, it is convenient to define a generating functional through introducing external sources jX(t) 
and jk{T) such that [21-23] 

Z\j*,j] = J D(a*a)e-'^*'=»'= / S)(a*a) exp{a^(/3)afc(/3) 
-IodT[al{T)dk{T)+K{a*a)-jl{T)ak{T)-al{T)jk{T)]} (42) 

where the sings and "— " in front of jj^^ refer to bosons and fermions respectively and Zo\j*,j] is defined by 

Zo[j*,j]= J D{a*a)e-°-'>''"' j D(a*a)e^("*'"'^'*'^') (43) 

in which 

I{a\a;j\j) = al((3)ak{l3) - dT[al{T)dk{T) 

+U!kal{T)akiT) -jl{T)ak{T) -4(T)jfc(T)] 

Obviously, the integral in Eq. (43) is of Gaussian-type. Therefore, it can be calculated by means of the stationary- 
phase method as will be shown in detail in Sect. 4. 

The exact partition functions can be obtained from the generating functional in Eq. (42) by setting the external 
sources to be zero 

^ = ^r,i] li-=i=o- (45) 

In particular, the generating functional is much useful to compute the finite-temperature Green functions. For 
simplicity, wc take the two-point Green function as an example to show this point. In many-body theory, the Green 
function usually is defined in the operator formalism by [21,29] 

Gki{TuT2) = ^Tr{e-P^T[au{T^)a+{T^)]} = Tr{eP^'^-^)T[au{T^)a+ {r^)]} (46) 

where < ri,T2 < /?, = — ^ InZ is the grand canonical potential, T denotes the "time" ordering operator, afc(ri) 

and af{T2) represent the annihilation and creation operators respectively. According to the procedure described in 
Eqs. (12)-(22). it is clear to see that when taking ti and T2 at two dividing points and applying the equations (5) 
and (6), the Green function may be expressed as a functional integral in the coherent-state representation as follows 

Gki{n,T2) = D{a*a)e-<<"' j ^{a*a)ak{Ti)a*i{T2)e'^''' '"^ (47) 

With the aid of the generating functional defined in Eq. (42), the above Green function may be represented as 

Gkl{Tl,T2) = X r ■ / X j*=j=0 (48) 

Z ^Jfc(n)dji(r2) 

where the sings " -|-" and " — " belong to bosons and fermions respectively. 
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III. THE COHERENT-STATE REPRESENTATION OF THERMAL QCD HAMILTONIAN AND ACTION 



To write out explicitly a path-integral expression of thermal QCD in the coherent-state representation, we first need 
to formulate the QCD in the coherent-state representation, namely, to give exact expressions of the QCD Hamiltonian 
and action in the coherent-state representation. For this purpose, we only need to work with the classical fields by 
using some skilful treatments. Let us start from the effective Lagrangian density of QCD which appears in the 
path- integral of the zero- temperature QCD [22, 25, 26] 

£ = i,{i^'^{d^ - igT'^Al) - m}V - - ^ {d^A^f - d^C^DfC^ (49) 

where T°- = X°-/2 is the color matrix, tp and tp represent the quark fields, are the vector potentials of gluon fields, 
C° and C" designate the ghost fields, 

F;, = d^A: - d.A; + gr'^A^Al (50) 

and 

Df = 5'^% - gr^'AI (51) 

For the sake of simplicity, we work in the Feynman gauge (a = 1). It is well-known that in this gauge, the results 
obtained from the above Lagrangian are equivalent to those derived from the following Lagrangian which is given by 
applying the Lorentz condition d'^A'^ = to the Lagrangian in Eq. (49), 

C = i,{ir{d^. - igT-A-) - m}^ - ^d^A-,d^^A-- - gt^-d^AlA^^'A-'^ 

_lg2j^abcfadej^biJ.j^cuj^dj^e _ dl^C°'d^C^ + g f^^d^CC^ A"^ ^ ' 

Here it is noted that the application of the Lorentz condition only changes the form of free part of the gluon Lagrangian, 

remaining the interaction part of the Lagrangian in Eq. (49) formally unchanged. The above Lagrangian is written 
in the Minkowski metric where the 7— matrix is defined as 70 = (3 and 7 = (3d [26]. In the following, it is convenient 
to represent the Lagrangian in the Euclidean metric with the imaginary time t = it where t is the real time. 

Since the path-integral in Eq. (42) is established in the first order (or say, Hamiltonian) formalism, to perform 
the path-integral quantization of thermal QCD in the coherent-state representation, we need to recast the above 
Lagrangian in the first order form. In doing this, it is necessary to introduce canonical conjugate momentum densities 
which are defined by [26, 30] 

nv = = *^7° = 
IT— — — n 

K = oi^ = -9tA; + gf^^A'^A^,, (53) 
= = -StC« + gr'^^c^A^o 

where the subscripts R and L mark the right and left-derivatives with respect to the real time respectively. With the 
above momentum densities, the Lagrangian in Eq. (52) can be represented as 

c = u^ dttp+ w^'dtA"^ + wdtC + dtC^u - n (54) 

where 

n = Ho + Hi (55) 

is the Hamiltonian density in which 

Wo = V;(7 • V + m) V + ^ (np2 _ V^^;! - WTf + (5" V^C" (56) 
is the free Hamiltonian density and 

Hi = igi^T-^^A-i, + gr'-{in-Al + diA^A^A^ - i^/V'^'V^'^^^^ (^7) 
X {AiAl - AIAI) + gf''^''{iWAl - diC^A^C^ ^ ' 
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is the interaction Haniiltonian density here the Latin letter i denotes the spatial index. The above Hamiltonian 
density is written in the Euclidean metric for later convenience. The matrix 7^ in this metric is defined by 74 = /3 
and 7 = —ipd [30]. It should be noted that the conjugate quantities 11" and n" for the ghost fields are respectively 
defined by the right-derivative and the left one as shown in Eq. (53) because only in this way one can get correct 
results. This unusual definition originates from the peculiar property of the ghost fields which are scalar fields, but 
subject to the commutation rule of Grassmann algebra. 

In order to derive an expression of the thermal QCD in the coherent-state representation, one should employ the 
Fourier transformations for the canonical variables of the QCD which are listed below. For the quark field [26, 30], 

i>{x,r) = J ^0j^y{^bMr)e'^^-'' + v^{m{P,r)e-'^"'] (58) 

r)=J [u' {p)K (P, r)e-^^^-^ + {p)d, {p, rje^^ (59) 

where u'^{p) and v^{p) are the spinor wave functions satisfying the normalization conditions u^~^{p)u^{p) = 
v'^^ {p)v'^ [p) = 1, bsip.T) and b*{p,T) are the eigenvalues of the quark annihilation and creation operators &s(p, t) and 
bf{p,T) which are defined in the Heisenberg picture, ds{p,T) and d*{p,T) are the corresponding ones for antiquarks. 
For the gluon field [26, 30], 



I _^e^^{k)[al{k,T)e'^-^ + al*{k,T)e-^^-^] (60) 



where e„( ) is the polarization vector and 



K{^^^)=i J (^y^^^WK(fc>T)e*'"^-«r(fc,^)e-*1 (61) 

which follows from the definition in Eq. (53) and is consistent with the Fourier representation of free fields. In 
the above, a^(fc, r) and a^*(fc, r) are the eigenvalues of the gluon annihilation and creation operators a^(fc, r) and 
a'^{k,T). For the ghost field, we have 

C'ix^r)=J 4i7i;7=[^«(.>)e^^^-><(.>)e----], (62) 

^"(-v) = / (^^[^»(^»^^'^"+^:(5»^-^'i' (63) 



n"(^'^) = ^ / (^V^t""^^"''^^"''"^"""^*"''^^""''"^' ^^^^ 



and 



^{^,r)=i J ^-^y^^[c„(g,r)e^f- -c:(9,r)e-^«-^ (65) 

where Ca{q,T) and c'^(q,T) are the eigenvalues of the ghost particle annihilation and creation operators Ca{q,T) and 
c^(g, r) and Ca{q,T) and c*(g, r) are the ones for antighost particles. 

For simplifying the expressions of the Hamiltonian and action of the thermal QCD, it is convenient to use abbre- 
viation notations. Define 

,er^ \ fbs{p,T), if e ^ +,\ 
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and furthermore, set a = (p, s, ^) and 



(27r)-3/2i;^(p), i/ 

E = E / (68) 

a se ^ 



Eqs. (58) and (59) may be represented as 

V(f,r) = E^'^o&a(T)e^^P^ 



Similarly, when we define 



a-(A:,r), if 9 = 



Afeik) = (27r)-3/2(2c.(^))-V2e;;(^), 
n^^,(fc) = i^(27r)-3/2[a,(g-)/2]V2eA(fc) 



Eqs. (60) and (61) can be written as 

Al{x,T) = j:A'^a„{T)e'6^-^ 



For the ghost fields, if we define 



a a9 

then, Eqs. (62)-(65) will be expressed as 

C"(f,r)=EGaC„(T)e^^^- 
C"^(f,T) =EG„4(r)e-^»^- 
^"(.x, t) =J2 n„c„(r)e*^^9-^ 
n"(£,r) = Enac;(r)e-^?- 

a 

Upon substituting Eqs. (69), (73) and (77) into Eqs. (56) and (57), it is not difficult to get 



(69) 



:a'i(k.T), if 9=+, , 



(71) 



and furthermore, set a = {k, c. A, 0) and 

E = E / d% (72) 



(73) 



G,(g-) = (27r)-3/2[2u;(g)]-i/2, 

n,(g-) = i^(27r)-3/2[a;(9)/2]V2, ^^^^ 

and furthermore set a = [q, a, 6) and 



(77) 
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Hoir) = J d^xHoix) = E^aea&*(r)&a(r) 
+ 5E'^a<(r)a„(r) +j:oJc,cl{T)ca{T) ^^^^ 

a a 

and 

Hi{t) = J (fxHiix) = E A{aPl)hl,{T)b0{T)a^{T) + E S(a/37)aa(T)a;3(r)a^(r) 

a/37 "/37 /'j-g-v 

+ E C'(a/37<5)aa(T)a/3(T)a^(T)a5(T) + E -C'(a/37)Ca(T)c;3(T)a^(T) ^ 

which are the QCD Hamiltonian given in the coherent state representation. In Eq. (78), the first, second and third 
terms are the free Hamiltonians for quarks, gluons and ghost particles respectively where 6^ = 9, Sa = {p^ +to^)^/^ is 

the quark energy, u>a = k is the energy for a gluon or a ghost particle. In Eq. (79), the first term is the interaction 

Hamiltonian between quarks and gluons, the second and third terms are the interaction Hamiltonian among gluons and 
the fourth term represents the interaction Hamiltonian between ghost particles and gluons. The coefficient functions 
in Eq. (79) are defined as follows: 

A{apj) = ig{2nf5\eo.Pa - 90/3 - e^k^)wl: (Pa)T«7^W^'f iP^K^ (^7). (80) 

(81) 



B{al3j) = ig{2TT)H\e^k^ + 0pkp + ^^fc^)r''^[n;^° (fc„) 
^^461^(^7) + ^ak"A1^g'^{ka)A^g^{kJ)]A^f^{k|3), 



and 



£)(a/37) = tg{2TTfS^9^q^ - 9^q^ - 0^fc^)/«''=GgJg;) 
X [n'epmAfe: ik,) - OckfGl^ (qfMte: (^7)]- 



(82) 



(83) 



It is emphasized that the expressions in Eqs. (78) and (79) are just the Hamiltonian of QCD appearing in the path- 
integral as shown in Eq. (42) where all the creation and annihilation operators in the Hamiltonian (which are written 
in a normal product) are replaced by their eigenvalues. 

To write the path-integral of thermal QCD, we need also an expression of action S given in the coherent state 
representation. This action can be obtained by using the Lagrangian density shown in Eq. (54). By partial integration 
and considering the following boundary conditions of the fields [20-22]: 



V'(:e,0) = xp{x), ■ip{x,0) = il){x), 
tp{x,P) = -tp{x), ip{x,l3) = -tp{x) 



A^(f,0)=^^(f,/3)=^^(f), 

n^(f,o) = n^(f,/?) = n^(f) 



and 



C"(f,0) = C''{x,p) = C"(x), C"(x,0) = C»(f,/3) = Cix), 
If (f,0) = ir(x,/3) = Ti (x), U"-{x,0) = U"-{x,f3) = n«(f), 

the action given by the Lagrangian density in Eq. (54) can be represented in the form 

S = lo dr J d3x{i[V.+ (x, T)V;(,x, r) - ijj+{x, T)V'(f , r)] 
+1 [n^(f , t)4 {x, t) - (f , r) (x, r)] 

+ ^[Ua{x,T)Ca{x,T) - Uaix,T)Ca{x, t) 

+ Ca {x, r)n„(f , r) - Ca{x, t) Tia {x, r)] - n{x, t)} 



(84) 



(85) 



(86) 



(87) 
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where the first relation in Eq. (53) has been used and the symbol " • " in ■0(a', t), ^^(.x, t) now denotes the 

derivatives of the fields with respect to the imaginary time t. It is stressed here that only the above expression is 
appropriate to use for deriving the coherent-state representation of the action by making use of the Fourier expansions 
written in Eqs. (58)-(65). On inserting Eqs. (58)-(65) into Eq. (87), it is not difficult to get 

S = -jyT{Jd^k{^[K{k,r%{k,T)-b*{k,r)b,{k,r)] + ^[d*{k,T)d,{k,T) 

-dt{k,T)ds{k,T)] + ^[al*{k,T)dl{k,T) - dl*{k,T)al{k,T)] + i[c:(fc,r) C„ (fc,T) (gg) 
{k,r)Ca{k,T)-cl{k,T)Ca{k,T)+cl{k,T)Ca{k,T)]} + H{T)} 

= —Se 

where H{t) is given by the sum of the Hamiltonians in Eqs. (78) and (79) and Se is the action defined in the 
Euclidean metric. It is noted that if one considers a grand canonical ensemble of QCD, the Hamiltonian in Eq. (88) 

should be replaced by K{t) defined in Eq. (2). Employing the abbreviation notation as denoted in Eqs. (66), (70) 
and (74) and letting stand for (oq,, ba, Ca), the action may be compactly represented as 

Se= f dT{Y,\K{r)oqo.{r)-ql{T)oq^{T)]+H{T)} (89) 
JO ^ 

where we have defined 

?a ° 9a = CLa- + b*J)a + OaC^Ca (90) 

It is emphasized that the 6'ct = ± is now contained in the subscript a. Therefore, each a may takes a+ and/or a~ as 
the first term in Eq. (90) does. 



IV. GENERATING FUNCTIONAL OF GREEN FUNCTIONS FOR THERMAL QCD 

With the action Se given in the preceding section, the quantization of the thermal QCD in the coherent-state 
representation is easily implemented by writing out its generating functional of thermal Green functions. According 
to the general formula shown in Eq. (42), the QCD generating functional can be formulated as 

ZL?-] =/i?(g*g)e-^*«/S)(g*q)exp{i[9*(/3)-g(/3) 

-q*{Q).q{Q)]-SE + S!^dT3*{T)-q{T)} ^ ^ 

where we have defined 

q* -9 = ^a* fla + Oababa + c* (92) 

and 

j*-q = + OaivlK + blrja + QCc, + 4Ca) (93) 

here ^a,ila and are the sources for gluons, quarks and ghost particles respectively and the repeated index implies 
summation. It is noted that the product q* ■ q defined above is different from the g* o q„ defined in Eq. (90) in the 
terms for quarks and ghost particles and the subscript a in Eqs. (92) and (93) is also defined by containing = ± 
. In what follows, we assign to represent the a with 6a = ±- According to this notation, the sources in Eq. (93) 
are specifically defined as follows: 

€ (94) 

c 

where the subscript a on the right hand side of each equality no longer contains 9a and the gluon term in Eq. (92) 
(l/2)a*Oa may be replaced by a„-a„+. The integration measures D{q*q) and D{q*q) are defined as in Eqs. (24) and 
(25). 

The generating functional in Eq. (91) is nonperturbative. Now we are interested in describing the perturbation 
method of calculating the QCD generating functional. Since the Hamiltonian can be split into two parts Ho{t) and 



Sa+ Sa) Sa~ 

Va+ = Va, Va- = 

Ca+ Cai Ca^ 
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Hj{t) as shown in Eqs. (78) and (79), the generating functional in Eq. (91) may be perturbatively represented in the 
form 

Z[j] = exp{- J^^ dr/i-,(^)}^°b1 (95) 

where is the generating functional for the free system and the exponential may be expanded in a Taylor series. 

In the above, the commutativity between Hj and Z'^lJ] has been considered. Obviously, the Z^[j] can be written as 

ZO[i] = Z°[e]ZO[r7]ZO[C] (96) 

where Zg[£,], Zg[r]] and Z^[(] are the generating functional contributed from the free Hamiltonians of gluons, quarks 
and ghost particles respectively. They are separately and specifically described below. 

In view of the expressions in Eqs. (91), (88) and (78), the generating functional Zg[l^] is of the form 



X jD{a*a) exp{/g(a^, ax;Cx,^x)} 



(97) 



where 



/g(a* , ax; Cx, Ca) = / rf^fciK(fc, P)ax{k, P) + a* (fc, 0)ax{k, 0)] 
-/(f dr/d3A;{i[a^(fc,r)aA(A^,r) -a^(fc,r)aA(fc,T)] +u;(fc)a^(^,r)aA(^,r) (98) 
-Cx{k,T)ax{k,T)-alik,T)^x{k,T)} 



and 



D{a*a)=Ulda*x{k)daxik), 

k\ 

= n ldal{k,T)dax{k,T). 

kXr 



(99) 



The subscript A in the above is now assigned to denote polarization and color. When we perform a partial integration, 
Eq. (98) becomes 



J9K,aA;a,eA)^= /d'fc<(fc,/3)«A(fc,/3) -/o'^dT/d3fc{a*(^,T)dA(fc, 
+uj{k)al{k, T)ax{k, t) - Cx{k,T)ax{k, t) - al{k,T)Cx{k, r)}. 

For the generating functional Z^ [r]] , we can write 

Z^irj] = J D{b*bd*d)exp{- J d^k[b%k)bsik) + d:ik)d,ik)]} 



X jT){b*bd*d) exp{7g(6*, b^, rf*, ds; v*s,Vs, V*s,Vs)} 



(100) 



(101) 



where 



Ig (6* , 6„ , 4 ; , , , r? J = / d^fci [6* (fc, p)b, (fc, p) + d: (fc, /3)4 {k, P) 

+b*{k,Ql)b,{k,Q)+d*{k,i))ds{k,Q)] - dT j (fk{^(k,T)b,(k,T) 

-b* (fc, T)bs (fc, r)] + i [dl {k, T)ds (fc, r) - (fc, t)4 (fc, r)] (102) 
+e{k)[b;{k,T)b,ik,T)+dt{k,T)d,{k,T)] - [r,;{k,T%ik,T) 
+bt{k,r)r],{k,T) +r]*{k,T)d,{k,T)+dl{k,T)v,k,r)]} 



and 



D{b*bd*d) = '[ldb*ik)dbsik)ddl(k)ddsik), 

ks 

D{b*bd*d) = n dbl{k,T)dbs{k,T)ddl{k,T)dds{k,T) 

ksT 



(103) 



in which the subscript s stands for spin, color and flavor. By a partial integration over r, Eq. (102) may be given a 
simpler expression 
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- lo / '^''HK (k, T)bs{k, t) + d; [k, T)d,{k, t) 

+e{k)[bl{k,T)b,{k,T)+d:ik,T)ds{k,T)] - [7?;(ft,T)6,(fc,T) 

+bt {k, r)r]s {k, r) + {k, T)ds {k, r) + rfj (fc, T)r],k, r)] }. 



(104) 



As for the generating functional .^c[C]) we have 



Z'eiC] = I Dilf-CCC*) exp{- / d^k[cl{k)-Ca{k) - <(_fc)Ca(fc)]} 

X jD{c*ccc*) exp{7c(c* , c„, c*, c„; C* , Ca, Cl, Ca)} 



where 



and 



/e(<, C,, C„; C, Ca,C, Ca) = / '^'fc5K(^./5)Sa(fc, - c:(fc, /3)Ca(fc, /3) 

+c:(fc,0)ca(^,0) -c:(fc,0)ca(fc,0)] -/o''dr/d3fc{i[c:(^,r) (fc,r) 

- cl {k, T)ca{k, r)] - i[<(fc, r)ca(fc, r) - c* (fc, r)ca(fc, r)] 
+a.-(fc)[c:(fc,T)ca(A^,r)_- <(! t)c„(^, r)] - [C* (fc, t)c„(^, r) 
+<(fc, T)Ca(fc, r) + Cik, r)ca{k, r) + c*(fc, T)C„(fc, r)]} 



£>(c*ccc*) = n(^c*(fc)cte„(fc)dc„(fc)dc*(A;), 

ka 

D{c*ccd*') = n dcl{k,T)dCa{k,T)dCa(k,T)dcl{k,T) 

kar 



(106) 



(107) 



in which the subscript a is a color index. After a partial integration, Eq. (106) is reduced to 

(C* , C„ , C* , Ca ; C , Ca , Cl , Ca) = / d^HK {k, /3)Ca (fc, /3) - C* (fc, /3) (fc, /?)] 

/(f dT/d3fc{c*(fc,T) C„^(fc,r) -C*(fc,T)c„(fc,r) +tj(fc)[c*(^,T)Ca(fc,T) ^^Qg^ 



-c:(fc,r)ca(fc,T)] - [C:(fc,r)ca(fc,T)_+<(fc,r)Ca(fc,T) +Ca(A;,T)ca(A;,r) 

+c:(fc,T)C(fc,r)]}. 

Here it is noted that all the terms related to the quantities c* and Ca are opposite in sign to the terms related to 
the c* and Co and, correspondingly, the definitions of the integration measures for these quantities, as shown in Eq. 
(107), are different from each other in the order of the differentials. 

The generating functional in Eqs. (97), (101) and (105) arc all of Gaussian- type, therefore, they can exactly 
be calculated by the stationary-phase method. First, we calculate the functional integral -^gi^]- According to the 
stationary- phase method, the functional Zg[^] can be represented in the form 

= I i?(a*a)exp{- / d^k[al{k)ax{k) + l'l{al,ayCx,W (109) 



where Ig{a\, a\; C^, $a) is given by the stationary condition 5Ig{a\, ay, CJ) Ca) = 0. By this condition and the boundary 
condition [20-22]: 

al{k,(3)=al{k), ax{k,0) = ax{k), (110) 
one may derive from Eq. (98) or (100) the following inhomogeneous equations of motion [23-26]: 

a\{k, t) + uj{k)ax{k, r) = ^x{k, r), 

al{k,T)-iv{k)al{k,T) = -^l{k,r). ^ ' 

In accordance with the general method of solving such a kind of equations, one may first solve the homogeneous linear 
equations as written in Eq. (33). Based on the solutions shown in Eq. (34) and the boundary condition denoted in 
Eq. (110), one may assume [23-26] 

ax{k, t) = [ax{k) + ux{k, T)]e-(fc>, 
alik, r) = [alik) + <(fc, T)]e-(^^)(-/') 
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where the unknown functions ux{k,T) and r) are required to satisfy the boundary conditions [23-26]: 

ux{k, 0) = uxik, p) = ulik, 0) = ul{k, p) = 0. (113) 
Inserting Eq. (112) into Eq. (Ill), we find 



Integrating these two equations and applying the boundary conditions in Eq. (113), one can get 

ux{k,T)=J^dT'e'^(^yUk,T'), 
Kik,T) = -/;dr'e-«(/3-')$*(fc,T'). 

Substitution of these solutions in Eq. (112) yields [23-26] 



ax{k,T) = aA(fc)e-W- + /;rfT'e— W(--')^^(fc,r'), 
al{k,T) = a^(fc)e'^(^)(^-'5) + /^^rfr'e'^(^)(^-^')^^(fc,T'). 



(114) 



(115) 



(116) 



When Eq. (116) is inserted into Eq. (98) or Eq. (100), one may obtain the (a^, a^; ^a) which leads to another 
expression of Eq. (109) like this 

Zg[^] = J D{a*a)e^p{- J d^k[al{k)axik){l - e^^^^^^) - a* (fc)e-'3'^(^) 

X dTe^Ck)r^^^k, t) - Jl^ dTe-^('^^^Cxik, r)axik)] (117) 
+ fo !^ dT2 J d^kCxik, n)e{n - T2)e-(^")(— =)a(fc, T2)}. 

When we set 

A = 1 - e-'^"^^), 

b = e-^'"^^) dTe^Ck)r^^^k, t), (118) 
/(a)=/,^dre-(^>ei(fc,r)aA(fc), 

by employing the formula denoted in Eq. (37), the integral over a^(fc) and ax{k) in Eq. (117) can easily be calculated. 
The result is 

Z0[e] = Z0exp{- Trfn [\t2 I d'k^^%k,n)Ali,{k,n-T2)^^:{k,T2)} (119) 
Jo Jo J 

where 

Z0= Y[[l- e-/3-(fe)]-i = [][!- e-^"^^^]-^ (120) 

k Xa k a 

is precisely the partition function contributed from the free gluons [21-24, 29] and 

A^^: (fc, n - T2) = gxx'S'^''' Ag{k, n - T2) (121) 

with 

Ag(fc, n - T2) = e{n - T2) - (1 - e/''^(^))-ie-'^(^'(^i-^=) (122) 

is the free gluon propagator given in the Feynman gauge and in the Minkowski metric {Note: in Euclidean metric, 
gxx' ^ ^<^AA')- III Eq. (119), the color index 'Vi" has boon explicitly written out and the A now merely designates 
the polarization index. In the other expressions, we still use A to mark the both of color and polarization indices for 
simplicity. When we interchange the integration variables n and T2 and make a transformation k — > —k in Eq. (119), 
by considering the relation 

Cxilr) = ^x{-k,T) (123) 
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which will be interpreted in the Appendix, one may find that the propagator in Eq. (122) can be represented in the 
form 

Ag{k,Ti - T2) = i[n6(fc)e-'^(^)l^i-^=l - n(,(fc)e'^(^)l^i-^=l] (124) 

where 

nt,{k) = (1 - e-^<''Y\ Mk) = (1 - e^^(^))-i (125) 

are just the boson distribution functions [20-23, 29]. 

Let us turn to the calculation of the functional integral in Eq. (101). Based on the stationary-phase method, we 
can write 

= J Dib*bd*d)exp{- J d^k[b:{k%{k) + d:{k)dM} (126) 
xexp{I°{bl,bs,dl,ds;vt,Vs,V*s,Vs)} 

where 1° (6*, 6^, dj, ds; yyj, ry^, Ty*,^^) will be obtained from Eq. (102) or Eq. (104) by the stationary condition 
SIq{b*,bs,d*,ds;r]*,r]s,r]*,r]g) = 0. Prom this condition and the boundary conditions [21-26]: 

btik,P) = -b*M, bs{k,0) = bM, ^27) 
dt{k,(3) = -d:{k), ds{k,0) = ds{k), ^ ' 



one may deduce from Eq. (102) or Eq. (104) the following equations [23-26]: 

(128) 



bs{k,T) +e(k)bs{k,T) = r]s{k,T), 
bl{k,T) - e{k)bl{k,T) = -t^HKt), 

ds{k, r) + e{k)ds{k, t) = 7]^(k,T), 
d*Ak,T) - e{k)dl{k,T) = -ril{k,T). 



Following the procedure described in Eqs. (111)-(116), the solutions to the above equations, which satisfies the 
boundary conditions in Eqs. (127) and the conditions like those in Eq. (113), can be found to be [23-26] 



(129) 



bs{k,T) = 6,(fc)e-^W- + /;rfr'e-^W(--')77,(ft,T'), 
b*{k,T) = -6:(fc)e^W(^-'5) +/^''dr'e^W(^-^')r;*(fc,T'), 

ds{k,T) = 4(fc)e-'^(^)^ +/„''dT'e-'^(^)(^-^')?7,(fc,r'), 
rf*(fc,r) = -d:(A^)e^(*^)(^-« +/^^dr'e'^(^)(^-^')77*(fc,T'). 

Substituting the above solutions into Eq. (102) or (104), we find 

I°{bt,bs,dl,d,,:rit.ns,frs,r]s) = /rf^A;{-e-/'-W[6:(fc)6,(fc) + d:(fc)d,(fc)] 

+ dTe-'^'^^^[r,*(k,T)b,{k) + r]*{k,T)ds{k)] ~ 6"^^^ dre^^^^ (130) 
X [6^ {k)r]s {k, t) + d; {k)r], {k, r)] } + B[r]; , % , 7?* , ??,] 

where 

B[v:,Vs,V:,Vs] ^U„^dT,J|^dT,Jd^k{ein-T,)e--Ck)ir^-r.) 

x[r]t{k,Ti)Tjs{k,T2) +rj:{k^n)r},{k,T2)] + e{T2_j- ri)e-^W(^=-^i) (131) 
X Ws {k, T2)vsik,n)+ fj* {k, T2 )r], {k, n )] } 

On inserting Eq. (130) into Eq. (126), we have 

Z°[,7] = A[r,:,r,„^:,?7je^[''^''-^*''?=1 (132) 

where 

AK, r^^r?:, TyJ = / Dib*b) cxp{- / d^k[b:{k%{m + e"^'^^^) 

+e-f'<^^bm jyre-C^^-Vsik^T) /^^^ rfre'^^^^)-,?: (fc, T)6.(fc)]} .^33) 
X / Did*d) exp{- / d^k[dt{k)ds{k){l + e-'^^C^)) + e-'^'^^'>dlik) 
X g dre<^)-%{k, r) - J,^ dre-^^^-fJUk, r)ds{k)]} 
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here the fact that the two integrals over {b* , b} and {d* , d} commute with each other has been noted. Obviously, each 
of the above integrals can easily be calculated by applying the integration formulas shown in Eq. (39). The result is 

A[vt,Vs,rj*s,Vs] = Z\ exp{-i dr^ J d^k{l + e^<^^)-' 

x{e--W(-i-2)[^*(fc,^^)^^(^^^2) +^:(^,ri)7?,(fc,r2)] (134) 

where 

Z0 = [][l + e-^^«]2 (135) 
k s 

which just is the partition function contributed from free quarks and antiquarks [21-23,29]. It is noted that the two 
terms in the exponent of Eq. (134) are equal to one another as seen from the interchange of the integration variables 
Ti and T2. After Eqs. (131) and 134) are substituted in Eq. (132), we get 

= exp{/(f dn /o^ dT2 J d^k{[0in - T2) - (1 + e'3-(fc~))-i]e-=(^)(-i--=) 
x[rjt{k,n)vs{k,T2) + rn{k,n)rl,{k,T2)] + [0{t2 - n) - (1 + e/^-W)"!] (136) 
xe-^W(-=--i)[r?;(fc, r2)?7.(fc, n) + r?:(fc, T2)^,(fc, Ti)]}}. 

When we interchange the variables n and T2 and set ^ ^ — fc in the second term of the above integrals and notice 
the relation 

V*s{k,T2)r],{k,n) = r]*{-k,Ti)r]s{-k,T2) (137) 
which will be proved in the Appendix , the functional [rj\ will eventually be represented as 

Zl[ri[ = Z° eMio dri J^^ dr2 J d^k[ri* [k, n) AJ«' (fc, n - T2)%. (fc, T2) .^3g. 
+r^Ak.T^)/^f{k,T^-T2)vAk,T2)]} 



where 



with 



A- (fc, n - T2) = 5'' Ag{k, n - T2) (139) 



A,(fc, n ~ T2) = ^[nfik)e-<^^\^'-^'\ - n/(fc)e"(^)l"i-"=l] (140) 
is the free quark (antiquark) propagator. In the above, 

nf{k) = (1 + e-^<''^)-\ nf{k) = (1 + e^^^^^)"^ (141) 

are the fermion distribution functions [21-23, 29]. 

Finally, 

let us calculate the generating functional .^cIC]- From the stationary condition (Co, Ca,c*,Ca; Ca) Co 5 Ca> Co) = 
and the boundary conditions: 

C:(fc,/?) =c:(fc), c:(fc,/3) = cUk), C„(fc,0) = Ca(fc), Ca{k,0) = Ca{k), (142) 

which are the same as those for scalar fields other than for the fermion fields [21], it is easy to derive from Eq. (106) 
or (108) the following equations of motion: 

Ca{k,T) +Uj{k)Ca{k,T) = -C,a{k,T), 

clik, t) - u;{k)c:{k, r) = C(^_^, r), 

Ca (k, t) + U){k)Ca{k, t) = Caik, t), 

K {k, t) - uj{k)cl{k, t) = -Cl(fc, t). 
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By the same procedure as stated in Eqs. (111)-(116), the solutions to the above equations can be found to be 

Ca(fc,r) = c,(fc)e-(^)- - /;dT'e-(^")(---')Ca(fc,T'), 

c„(fc,r) = c„(fc)e-W-+/;dT'e-W(---')C_Jfc,r'), ^ ' 

Upon substituting the above solutions into Eq. (106) or Eq. (108), we find 

I°{cl,Ca,cl,Ca;Ca,Ca,Cl(a) = / d^fcie-^^-^^^) ^ (fc)c, (fc) - C*(fc)c„(fc)] 
+e-'3-('=) J^' dT6-('=)-[c*(fc)Ca(fc, t)_+ cmCaik, t)] 

- dre--Ck^^[Qik,T)ca{k) +Cik,r)cam (145) 
- lo dn !^ dT2 J d^k9{n - r2)e--(^)(— [C(^, n)Uk, r^) 
-C{k,ri)Ca{k,T2)]. 

On inserting the above expression into the following integral given by the stationary-phase method: 

ZclQ = J D{C*CCC*) exp{- J d^k[c:{k)Ca{k) - cl{k)Ca{k)] + I°{cl, Ca,C*a,Ca; Ca,Ca,ClCa)}, (146) 

and applying the integration formulas in Eq. (39) , one can get 

ZO[C] = ^° exp{/o'^ dn I^' dr,J d'k[e{n - T,) - (1 - e/5-(S))-i]e-W(-i-^) 



where 



Z° == [|[1 -e-'^'^(^)]2 (148) 

k a 



is just the partition function arising from the free ghost particles which plays the role of cancelling out the unphysical 
contribution contained in Eq. (120). If we change the integration variables in Eq. (147) and considering the relations 

C(fc,T) = -U-k,r), Ca{k,T) = -ta{-k,r) (149) 
which will be interpreted in the Appendix, Eq. (147) may be recast in the form 

ZO[C] = ^0 exp{/,'^ dn !^ dT2 J d^kl^lik, n)Af {k, n - r2)C' 0, r^) .^^q. 
-C(fc,ri)Ar'(fc,n-T2)Ca'(^,r2)]} 



where 

Ar'(fc,Tl - T2) = 5'"''Ag{k,T^ - T2) (151) 



here Ag(fc, n — T2) was written in Eq. (124). 

Up to the present, the perturbativc expansion of the thermal QCD generating functional in the coherent-state 
representation has exactly been given by the combination of Eqs. (95), (96), (119), (138) and (150). In the derivation 
of the perturbation expansion, as one has seen, to obtain the final expressions of the propagators shown in Eqs. (124), 
(140) and (151), it is necessary to use the functional properties and relations for the external sources as denoted in 
Eqs. (123), (137) and (149). As a result of the derivation of the generating functional, the partition function for 
the free system has simultaneously been given by the combination of Eqs. (120), (135) and (148). The partition 
function for the interacting system can be calculated in the way as shown in Eq. (45). Here it should be noted that 

the differential d/dj{T) in Eq. (95) represent the collection of the differentials i5/(5^^*(/c, r), 6/6^'^{k,T), —d/Sr]s{k,T), 
S/Sr]*ik,T), -S/Sr],ik,T), d/5r}*ik,T), -5/5C,a{k,T), 6/5Q{k,T), -5/5C,^{k,T) and S/dfaik-.r). Ordinarily, the gen- 
erating functional of thermal QCD represented in the position space is used in the literature. In the Appendix, it 
will be shown that this generating functional can readily be derived from the generating functional described in this 
section. 
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V. RELATIVISTIC EQUATION FOR QQ BOUND STATES 



With the generating functional given in the preceding section, we are ready to derive the relativistic equation 
for qq bound states at finite temperature. It is well-known that a bound state exists in the space-like Minkowski 
space in which there always is an equal-time Lorentz frame. Since in the equal-time frame, the relativistic equation is 
reduced to a three-dimensional one without loss of any rigorism, in this section we only pay our attention to the three- 
dimensional equation which may be derived from the equations of motion satisfied by the following two-" time" 
(temperature) four-point Green function [15-17] 

e(a/3;7'^;Ti -T2) = Tr{e'='(f^-^)T{iV[6„(Ti)&+(Ti)]iV[6^(T2)6+(r2)]}} 

;r{iV[6„(n)S+(ri)]7V[6,(r2)6+(r2)]})^ ^^^^^ 

where the symbol (• • •)^ represents the statistical average and N symbolizes the normal product whose definition can 
be given from the corresponding definition at zero-temperature by replacing the vacuum average with the statistical 
average [17] 

^^[?a(ri)S^(T2)] = r[6„(Ti)6+(r2)] - 5„^(ti - T2) (153) 

where 

5«/3(ti - r2) = (t[6„(ti)6^(t2)])^ (154) 

is the quark or antiquark thermal propagator. The normal product in Eq. (152) plays a role of excluding the 
contraction between the quark and the antiquark operators from the Green function when the quark and antiquark 
arc of the same flavor. Physically, this avoids the qq annihilation that would break stability of a bound state. 
Substituting Eq. (153) in Eq. (152), we have 

g(a/J; 7(5; n - rs) = G(a/3; 7^; n - T2) - So^nS^^ (155) 

where 

G(a/3; 7^; n - T2) = (T{6„(Ti)6+(ri)6^(r2)6+(T2)})^ (156) 

is the ordinary Green function and, Sa^ and S-y^ are the equal-time quark (antiquark) propagators. Obviously, in 
order to derive the equation of motion satisfied by the Green function ^(a/3; 7i5; n — t-i), we need first to derive the 
equation of motion for the Green fimction G{af3; jS; ti — T2). 

Let us start with the generating functional in Eq. (91). As shown in Sect. 4, by partial integration of the second 
term on the right hand side of Eq. (89), the generating functional may be written in the form 

m = J D{q*q)e-'^'-^ J D{q*q) exp{q* (/3) • q{P) 

-SE + I,^dTf{r)-q{r)} ^ > 

where 

Se = / dT{J2 Ql{r) o g„(T) + H{t)} (158) 

here H{t) was given in Eqs. (78) and (79). First, we derive an equation of motion describing the variation of the qq" 
four-point Green function with the "time" variable t\. For this purpose, let us differentiate the generating functional 
in Eq. (157) with respect to 6* (ti). Considering that the generating functional is independent of V^{t\) and noticing 
the expressions given in Eqs. (158), (78), (79) and (93), one may obtain 

= / D{q*q)e-i^ '> j ^{q* q)[-b,M) - eM^{r,) 
-Yl Aapxbp{Ti)a\(Ti) + 0aria{Ti)]cxp{q* {P) ■q{P) - Se (I59) 

-J,^dTr{r).q{T)} = 0. 
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When the ^^(ti) and a\{Ti) in the above are replaced by the functional derivatives 9aS/5ri'^{Ti) and ^/5j^(ri) 
respectively and multiplying the both sides of Eq. (159) with 9a, the above equation can be written as 



d S „ 5 (^2 

pX 

Then, we differentiate the above equation with respect to the sources r//3((ri), giving 



, + OaSa . , + VMp^(apA) ^ - r]a{Ti)}Z[j] = 0. (160) 



f A (161) 



+^a/3-r?„(Ti)3j-.}Z[i]=0. 



Furthermore, successive differentiations of Eq. (161) with respect sources r]*{T2) and ris{T2) yield 

U d S \ I 5j 

IvdTi <5j)* (ri) ' (5r;^(Ti)(5r;*(r2)5775(T2) " " i5'7S(''"i)i5'73(''"i)i5'7^(''"2)i5'75(t-2) 

+ ^^«^^"^("^'^)*'7;(i-i)*w(ri)5r,*(r2)5r,5(r2)5j5;(ri) + ^"Z' i77;(r2)6r,ii(r2) (162) 

Similarly, when differentiating Eq. (157) with respect &/3(ri), one may obtain 

d 5 „ S - - " 

Subsequently, On differentiating the above equation with respect to '7^(ti), we get 



{j- , , - eaEB , , - Y^aOaAiafiX) ^ , , - r/S(ri)}Z[j] = 0, (163) 



{s5^(3?r3^;W^ ~ ^'^^'^SrJrf)Sr,0(Ti) Sr,'Jri)Sr,t{Ti)Sj*(ri) 

-<5a;3 + r?^(n)3^Vb1=0. 
Finally, successive differentiations of the above equation with respect to the sources r?*(T2) and r]s{T2) give rise to 

<5';S(ti) ^ li-ri (5j);3(ti) ' Sr]!,{'^2)Sr]s{T2) P P Sri^{Ti)Sr]fj{Ti)Sri*{T2)Sr]s{T2) 

P 0peaA{al3X) g^^, (ri)57;„(ri)<5r,;(T2)5%(T2)5,jJ(Ti) ^ '^"^^ <5r,; (rs),?^)^ (ra) (165) 



A, 

t)' 



+(5^^(5(ti - 7-2)5^.(^;)i-^^(^^;, + »?^(ti) 5^J(ri)54(T2)d-%(r2) i^I^'l " 0- 

Adding Eq. (161) to Eq. (164), then multiplying the both sides of the equation thus obtained with —6adj3 and 
finally setting the external sources rja — VP = 0) but remaining the gluon source j\ ^ 0, we get 

(;^ + OaSa - e0ep)Si^^ + Y,[A{ap\)5p, - A{al3\)5ap\j^^S% = (166) 

where 

= -^(^'^(^Psr,*Jn)6i{n) ^<=''^=° ^^^^^ 
is the quark (antiquark) equal- time propagator in the presence of source j\. If we define 

H{a0;pa;Tiy^ = {-j- + ^a£a - 6pe0)5ap5p^ + ^/(a/J;paA)-— - (168) 

where 

/(a/J; paX) = A{apX)dpa - A{af3X)dap, (169) 
Eq. (166) can simply be represented as 
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Y,H{ap;pa;ny^Sj,^, = 0. (170) 

per 

When summing up the both equations in Eqs. (162) and (165), then multiplying the equation thus obtained with 
OaO^O-yOs and finally sotting all the sources but j\ to be zero, one may get 

pa\ ^' (171) 



where 



and 



(5r?5(Ti)(5ry^(ri)(5ry*(r2)(5%(r2) 



are respectively the two-" time" four-point thermal Green function and the quark or antiquark thermal propagator 
in presence of source j\. When the source j\ is turned off, Eqs. (172) and (173) will respectively go over to the Green 
fmiction in Eq. (156) and the propagator in Eq. (154). It is noted that due to the restriction of the delta function, 
the propagators in Eq. (171) are actually "time" -independent. With the definition in Eq. (168), Eq. (171) may be 
represented as 

^H{ap;pa;Tiy'G{pa;-r5;n - ra)^'^ = -6{n - T2)S{ap;-fSy' (174) 

pa- 

where 

S{aP; ^Sy^ = (5„55^^^ - S^^Si^^^- (175) 

Acting on the both sides of Eq. (155) with the operator H{a(3; pa;Tiy^ and using the equations in Eqs. (170) and 
(174), we find 

Y,H{aP;pa;Tiy^g{pa;-fS;Ti - T2y^ = J2 H{aP; pa;Tiy^G{pa;'y6;Ti - ra)-''^ 

per (176) 

= -S{n-T2)S{a(3;jSy^ 

This indicates that the equation of motion satisfied by the Green function Q{a(}; j5; t\ — T2) formally is the same as 
the one shown in Eq. (171). Therefore, in the case that the source jx vanishes, we can write 

(^ + ^a£a ~ 9i3£f3)g{a(3] 7(5; Ti - T2) 

= -^(ri - T2)5(a/3; 7,5) - E P^A)e(paA; 7^; n - ra) (177) 

where 

(5 



g{pa\; 7i5; n - T2) = jj^j^Q{p(J] 7(5; n - T2)-'^ 
T{7V[6,(Ti)6+(n)k(n)]A^K(r2)6+(r2)]}; 



(178) 



/3 

and 



5(a/3; 7^) = S^^S^p - S^^S^s = -(M^, M|]-)/3. (179) 

It is noted here that similar to the definition in Eq. (153), the normal product N[bp{Ti)b^ {Ti)ax{Ti)] in Eq. (178) is 
defined as 

7V[6,(Ti)6+(ri)aA(Ti)] = T^pin)b+{n)ax{ri)] - A{pctX) (180) 
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where 

K{pa\) = (T[6,(ri)6+(Ti)aA(Ti)])^ . (181) 

Substituting Eqs. (153) and (180) into Eq. (178), wc have 

g{p(j\; -i5; n - T2) = G{paX; jS; n - ra) - A{pc7X)S^s (182) 

where 

G{paX; 7^; n - T2) = (T{bp{n)bt{n)ax{ri)bj{r2)bj{T2)}) (183) 



is the ordinary five-point thermal Green function including a gluon operator in it. 

By the argument as mentioned in the Appendix (see Eq. (A9) or (A18)), it is easy to prove that the Green functions 
g{aP; 7(5; n — T2) and G{pcrX; 7^; n — T2) are periodic. Therefore, we have the following Fourier expansions: 

g{ap; j6; r) = i E G{a/3; j6; a;„)e-'-"^ 

QipaX; 7^; r) = i E g{paX; jS; a;„)e-^-''- ^^^^^ 

n 

where t = ti — T2 and w„ = ^5^. Upon inserting Eq. (184) into Eq. (177) and performing the integration 
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5 1-0 dre^'^"'^, we arrive at 



{iu)n - OaSa + 9/3ei3)g{aP; 7(5; w„) 
= -S{ap; 7(5) + E /M; P^A)g(paA; 7<5; tOn). (185) 

It is well-known that the Green function ^(pcrA; 7(5; a;„) is B-S (two-particle) reducible [15-17]. Therefore, we can 
write 

/(a/3; paX)g{p(jX\ 7(5; a;„) = ^ X(a/3; /iz^; uJn)g{pu; -f5; a;„) (186) 

Arp fiiy 

where K{a(3; ni^;u)n) is called interaction kernel. Thus, Eq. (185) can be written in a closed form 

{iu)n - Oa£c, + O0S0)g{a(3; 7(5; a;„) = -S{a(3; 7(5) + ^ ii'(Q!/3; pu; ujn)g{p.v\ 7^; (^„). (187) 

Now, let us turn to the equation satisfied by qq bound states. This equation can be derived from Eq. (187) with 
the aid of the following Lehmann representation of the four-point Green function which may be derived by expanding 
the time-ordered product in Eq. (152) and then inserting the complete set of qq bound states into Eq. (152) [15-17, 
31], 



(7'^)Xmn(<2/3), , , 

g(a(3;-/d;uji) = -el" 2^A„„| — ■ , , ^ 1 (188) 

where 

Xnmiap) = (m 7V[M+] n) (189) 



which is the transition amplitude from the state with energy En to the state with energy E^n and 

A„„ = e-'='^"-e-^^'". (190) 

Upon substituting Eq. (188) into Eq. (187) and then taking the limit: limit^,,^£;„^(ia;; — Enm), we get the following 
equation satisfied by the transition amplitude 

{Enm - Oa£a + O()£f3)Xnm{a0) = ^ K{q.(3] 7(5; Enm)Xnm{l5) (191) 

7(5 
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where the fact that the function S{aj3] ^5) has no bound state poles has been considered. If we take | m) to be the 
vacuum state | 0) and set E = Eno and Xn{o:l3) = (o N[bab'^]^ nj, we can write from the above equation that 

{E - + eps,3)Xniap) = Y,K{a(i;^5;E)xn{l5). (192) 

7(5 

where the subscript n in En has been suppressed. This just is the equation satisfied by the qq bound states at finite 
temperature. 

Since the index a contains 6a = ±, Eq. (192) actually is a set of coupled equations for the amplitudes Xn{oi'^P~), 
X„(q^~/3^), Xn{<^^ (3^) said Xn(Q!~/3~). Following the procedure described in Rcfs. (16) and (17), one may reduce 
the above equation to an equivalent equation satisfied by the amplitude of positive energy. We do not repeat the 
derivation here. We only show the result as follows: 

[E - e{ka) - e{k0)]iP{aP; -E) = ^ V{aP; jS; E)^{^5; E). (193) 

where il^{a(3\ E) = Xn(ct"^/3~) and V{af3', 7(5; E) is the interaction Hamiltonian which can be expressed as 

y(a/3;75;£) = ^y(")(a/3;7'5;i?), (194) 



71=0 



in which 



l/(°)(a/3;7<5;£;) = K++++{aP;^5;E), (195) 

yW(a/3;7<5;£;)= V y K++'^>^i''fi-.P^-.E)K^,++{pa-n5-,E) 

E-ae{k,)-he{K) 



V^'^\al3;^5;E) 



{E-ae(k,)-he(k,))(E-ce(k,)-de(k.)) ' (197) 



here a, 6 = ±, and 



K++++{a(3;^5;E)=K{a+l3-;^+5-;E), 

K (a/3;7<5;£:) = K {a- (3+ ; 5+ ; E) 

K+_+_{ap; 7(5; E) = K {a+ P+ n+ 6+ ; E) , 
K_+_+{ap; 7(5; E) = K{a-p- ; ^-5-;E). 



(198) 



VI. CLOSED EXPRESSION OF THE INTERACTION KERNEL IN THE EQUATION FOR QQ BOUND 

STATES 

In this section, we are devoted to deriving a closed expression of the interaction kernel appearing in Eq. (192) 
and defined in Eq. (186). For this derivation, we need equations of motion which describe evolution of the Green 

functions Q{af3; "/6; ti — t^) and ^(q:/3(t; 7(5: ri — t^) with time T2. Taking the derivatives of the generating functional 
in Eq. (157) with respect to 6^(t2) and 65(t2) respectively, by the same procedure as described in the derivation of 
Eq. (160), one may obtain 

, ^, , + g-yg-y , ^, , 6^QpA{-ip\)——A:——— - r]^(Ti)\Z\i\ = 0. (199) 

and 
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^T-T^ - ^^^^^^ -^M.^MA) f - r7l(r2)}^L7l = 0. (200) 

Performing differentiations of Eqs. (199) and (200) with respect to the sources r]s{T2) and r]*{T2) respectively, we get 

and 



Furthermore, by successively differentiating Eqs. (201) and (202) with respect to the sources r?a(ri) and r]j3{Ti), one 
obtains 

r 5^ / d ^ \ ^, , \ e £ — 

+ S O-fdxAijpX) 5^* (n )<5»7a (ri )5r;; (T2)5r;5 (r2 {t2 ) + ^T* ^^tSI^^FwI^ (2^^) 



Act 



■/A 

-(5/3^(5(ri - 7-2) j^,(^^'^j^^(^^) - ^7 (^2 ) g^, (ri (ri )6r;, (r2 ) }^ ^1 = 



and 



5';J(Ti)(5»)/3(Ti)<5r/*(T2) Sris{T2)' ° ° Sri^{Ti)Srii3{Ti)Sri*{T2)Sris{T2) 



- J2 0s9^A{aSX) SrAri)Svn(ri)Sv:,(T2)Sv.{r2)Sjl(r2) ~ Sr,'Jr2)Sr,0{r2) (204) 
+'5a5(5(Ti - T-2)3^T^7;j5;jJ(77) + V*5 (^2) rf,,;(ri)5,,^(Ti)5r,;(r2) }^b1 = 0- 



Let us sum up Eqs. (201) and (202) at first, then multiply the both sides of the equation thus obtained with —0-y9s 
and finally set all the sources but the source j\ to vanish. By these operations, we get 

Y,H{-fS;pa;T2y'Sj,', = (205) 

ptT 

where 

¥(7<5; pa; ra)^'^ = (31- + ^7^7 " 9s£6)SjpSsa - V f{paX; iS)——- (206) 
dr2 ^ ojy^{T2) 

in which 

/(paA; 7(5) = A(aa),5^p - A(7pA))(55ct = -/(7<5; paA) (207) 

and S^p^ was defined in Eq. (167). 

When we sum up Eqs. (203) and (204), then multiply the both sides of the equation thus obtained with OaO^OjOs 
and finally set all the sources but the source j\ to be zero, according to the definitions in Eqs. (172) and (173), it is 
found that 

Y,H{-fd;pa;T2y'G{al3;j5;Ti - rs)^'^ = <5(ti - T2)[Sa5S^(}{T2 - nY' - S^jSasin - ra)^^]. (208) 

pa- 
in order to derive the equation of motion satisfied by the Green function G{XTa; jd; ti — T2) defined in Eq. (183), 
we may take the derivative of Eq. (208) with respect to Jxiri). In this way, we get 

J2h{i5; pa;T2y'G{apX;^5;Ti - ra)^'^ = 5{n - T2)[5asA{lPp;T2 - nY' - 6p^K{a6p;Ti - rs)^'^] (209) 

pa 

where 
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A(7/3p; T2 - T^y^ = — — S^^(t2 - T^y^ , (210) 
A(a,5p; n - ra)^^ = -jr^Sasin - ra)^'^ (211) 



and 



G(a/?A;7(5;Ti - T2) = — — G(a/3; 7^; n - ra)^'^ (212) 

Acting on Eqs. (155) and (182) with the operator H{'y6; pa;T2) and employing Eq. (205), we find 

J2h{iS; pa; T2y^g{a(3; pa; n - ts)^'^ = ^^^(7'^; P^; T2y^G{a(3; pa; n - rs)^'^ (213) 

pa pa 

and 

Y,H{j5;pa;T2y^g{aPX;pa;n - ra)-''^ = Y^^i^S; pa;T2y^G{a(3X; pa;n - T2)^^ (214) 

/9(T pa 

The above two equahties further indicate that the equations of motion satisfied by the Green functions G{a(3; pa; t\ — 
T2y^ and Q{af3\; pa;T\ — T2)^^ are formally the same as those for the ordinary Green functions G{a(3; pa;Ti — r2)-'^ 
and G{aPX; pa;Ti - T2)-'^ respectively. Upon inserting Eqs. (208) into Eq. (213) and Eq. (209) into Eq. (214) and 
turning off the source jx, noticing the definition in Eq. (206), we derive the following equations 

(3^ + ^7^7 - Oses)g{aP;-fS;Ti - T2) = S{ti - T2)[SasSjp{T2 - n) 

-Sfs^Sc^sin - T2)] + E g{aP;XTa;n - T2)f{XTa;-fS) (215) 

\ra 



and 



+ ^7^7 - 05£5)G{(^(3p; lS; Ti - T2) = S{ti - T2)[5a5H'lf3p; T2 - Ti) 

-60-yA{aSp; n - T2)] + E 0{af3p; Xra; n - T2)/(Ara; 7^) (216) 

\t(t 



where some indices have been changed for convenience, 

A(7/3p;t2 - ri) = {T[%{T2)b+{n)ap{n)]) ^ 
A{a6p;Ti - T2) = (r[6c«(Ti)6|(r2)ap(ri 



(217) 

which are given by Eqs. (210) and (211) with setting j\ = and 

Q{Xtp; 7(5cr; n - T2) = ^^-.^^^^g^-^^^^) ^(At; 7(5; n - ra)^^ 

= (T{iV[5;,(Ti)5+(n)i(ri)]iV[6,(r2)6+(T2)a<,(r2)]})^ ^^^^^ 

is the six-point Green function including two gluon operators in it. According to the definition in Eq. (180), we have 

Q{Xtp; jSa; n - T2) = G{Xtp; jSa; n - T2) - A(ATp)A(7fo) (219) 

where 

G{XTp;-f5a;n - T2) = (r[&A(ri)6+(Ti)ap(Ti)6^(T2)6+(T2)a,(T2)])/3 (220) 

is the ordinary six-point Green function. It should be noted that due to the restriction of the delta function, the 
terms in the brackets on the right hand sides of Eqs. (215) and (216) actually are "time" -independent. 

It is easy to see that the Green functions Q{a(3;XTa;Ti — T2) and Q{XTp;^5a;Ti — T2), as the Green functions 
Q{af3; ^5; t\ —T2) and Q{afip; 7(5; n —T2), are periodic. Therefore, by the Fourier transformation, i.e. by the integration 
dre"^""^, noticing d/dT2 = —d/dr, Eqs. (215) and (216) will be transformed to 
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+ E Ara; a;„)/(ATa; 7(5) (221) 

At(T 



where S{ap; jS) was defined in Eq. (179) and 



{itOn + O-ySj - 9ses)Q{al3p\ 7(5; u)„) = R{a(3p\ ^5) 

+ E Q{a(ip\ Ara; a;„)/(Ara; 7<5) (222) 

\t(t 



wliere 



R{al3p\ 7^) = SccsKllip) - Si3jA{aSp) 

wliich is "time" -independent. 

Now we are ready to derive the interaction kernel. Acting on the both sides of Eq. (186) with {ioJn + d-ySj — ds^s) 
and using Eqs. (221) and (222), one gets 

j:Kiad: fii^:LOn)S{piy:-fd) = ^ f{aP; Xtp)R{\tp;-,S) + E E /("/?; Arp) 
x0(Arp; ^77ct; w„)/(^r?CT; 7(5) - E E -^(q;/?; MZ^; w„)^(a/3; ^r?c7; a;„)/(AT(T; 7^). 

Operating on the both sides of Eq. (186) with the inverse of Q{p,v;^5]UJn)i we have 

K{al3; 7(5; ^n) =Y.Y. -^("^S Arp)e(Arp; IJ,u;Wn)Q ^(/xz/;7^;a;„). (225) 

7(5 Arcr 

Upon substituting Eq. (225) onto the right hand side of Eq. (224) and acting on Eq. (224) with the inverse 
S'~^(/iz/;7(5), we eventually arrive at 

Kiaf3; jS; i^) = E{ E /("/5; \rp)R{\Tp; H + E E f{<^P\ XTp)g{XTp; ^w, E)fi^w, H 

- E E E E /(«/?; Arp)e(Arp; «^; E)g-^ («c; nO; E)g{ne; ^rja; E)f{^7^a; H}^"' (m^^; 7-5) ^ ' 

where cj„ has been replaced by E. This just is the wanted closed expression of the interaction kernel appearing in 
Eq. (192). In accordance with Eq. (186), the last term in Eq. (226) can be written in the form 

E E E ^("/^; P''' E)QiP^-^ ^V, E)K{^m E)S-\iJiv- 75) (227) 

pa ^ri pv 

which exhibits a typical B-S reducible structure [17]. Therefore, the last term in Eq. (226) plays the role of cancelling 
the B-S reducible part contained in the other terms in Eq. (226) to make the kernel to be B-S irreducible. If we use 
the above expression in place of the last term in Eq. (226) and acting on Eq. (226) with S{'yS; p,v), we obtain from 
Eq. (226) an integral equation satisfied by the kernel K{aj3; 7(5; E). Define 

7^(a/3; 7(5) = ^ /(a/3; \tp)R{\tp; jS) (228) 

Xrp 

and 

e(a/3; 7(5) = ^ ^ /(a/3; XTp)g{XTp; ^rja; E)f{Cna; jS), (229) 

Xrp ^rj(T 

the integral equation can be written in the matrix form as follows 

KS = n+Q- KgK. (230) 

For comparison with the kernel in Eq. (226) and for convenience of nonpcrturbativc investigations, we would like 
to show the corresponding closed expression given in the position space without giving derivation. This kernel can be 
obtained from the kernel in Eq. (226) by making use of the inverse of the Fourier transformations written in sect. 3 
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or derived from the generating functional represented in the position space (see Appendix ) by completely following 
the procedure as described in this section. The kernel is represented as follows: 

K{xiX2;yi,y2]E) = J d^zid^Z2{Tl{xiX2;zi,Z2) 

+Q{xi,X2]Zi,Z2]E) -V{xi,X2]Zi,Z2]E)}S~'^{zi,Z2-,yi,y2) ^ ' 

where TZ{xiX2; ^i, ^2), Q{xi,X2\z, , Z2; E) and V{xi,x2;zi, Z2; E) are separately described below. 
The function 'R,{x\X2;z\,Z2) can be represented as 

2 

Tl{xx, X2;zuZ2) = J2 ^7'^'li^'"{xi,X2;zuZ2) (232) 

i=l 

in which 

fir = igiHT,^, ^2 = igiHTl (233) 
with Tf = X"/2 and = —\°'*/2 being the quark and antiquark color matrices respectively and 

7e«"(f 1, fa; 21, 5*2) = (53(fi - zi)^lKl^{xi I X2, Z2) + 5^{X2 - ^'2)72 A^(^i I ^i, ^1) (234) 
here ^'^{xi \ xi,yi) and A.^{xi \ X2,y2) are defined as 



(235) 



Alix^ I xi.yi) = {T[A^^ix,,n)iP{xi,Ti)-iP{yi,Ti)])(3, 

A™(X, I X2,y2) = (T[A«(x„Ti)V^'=(.f2,Ti)/(y2,Ti)])/3 

which are time-independent due to the translation-invariance property of the Green functions. 
The function is of the form 

2 

Q{xuX2; zi,Z2; E)=Y^ ^T^fM, h \ fi, fa; ^i, ^2; E)^- (236) 

in which 

nr = igiHTt. nr = igi^l\%, (237) 
Q'^^{xi, Zj I x\,X2;zi,Z2\E) is the Fourier transform of the Green function defined by 

Glt{^i,Zj\xi,X2\ZuZ2;Ti-T2) _ _ 

= {T{N[A-{xi,n)i,{xuri)rix2, ri)]Ar[A^(£,-, r2)^(fi, r2)/(f2, r2)]})0 
The function V{xi, X2; Zi, Z2; E) is expressed by 

T>{xi,X2;zi,Z2;E) = / H d^Ukd^Vk ^TG}^ i^i I xi,X2;ui,U2;E) 

k=l i,j=l (239) 

xg-'^{ui,U2;vi,V2; E)gi^^''{zj \ vi,v2; zi, Z2; E)n^'' 

in which ^^*^"(fi | f 1, fa; wi, W2; £;) and e^''^'(fj I •yi, V2; 01, ^2; -E) are the Fourier transforms of the following Green 
functions 

^^^'^"(fi I fi,f2;Mi,«2;Ti - T2) _ /240) 

= (T{Ar[A«(fi,Ti)V;(fi,ri)r(f2,ri)]7V[V^(5i,T2)V'(52,r2)]});3 ^ ' 



(238) 



and 

Qi^^^{Zj I Vi,V2]Zi,Z2]Ti -T2) 



{T{N[i, {V, , Ti)r iV2 ,n)]N [Ai iZj,T2Mz,,T2)^ {Z2 , T2)]}} 



(241) 



The S ^{zi,Z2;yi,y2) in Eq. (231) is the inverse of the function defined by 

S{xi,X2;Zi, Z2) = (5^(fi - Zi)-ilSp{x2 - Z2) + S^{X2 - Z2)j2SFixi - Z 1 ) (242) 

in which Sf{xi — Zi) and Sp{x2 — Z2) are the equal-time quark and antiquark thermal propagators respectively. It is 
clear that there is one-to-one correspondence between the both kernels written in Eqs. (226) and (231). It is noted 
that the interaction kernel derived in this section is nonperturbativc because the Green functions included in the 
kernel are defined in the Heisenberg picture. Perturbative calculations of the kernel can easily be done by using the 
familiar perturbative expansions of Green functions as illustrated in the next section. 
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VII. ONE GLUON EXCHANGE KERNEL AND HAMILTONIAN 



In this section, we would like to show the one-gluon exchange kernel given by the expression in Eq. (226). In 

the lowest order approximation of perturbation, only the first term of the scries in Eq. (194), i.e., the kernel 
K{a^(3~;^^d~;E) represented in Eq. (195) is needed to be taken into account in Eqs. (192) and (193). In the 
lowest order approximation, as seen from Eq. (223), the first term in Eq. (226) vanishes because the function 
R{Xt p: "f~^ S~ ) gives no contribution to the kernel owing to the expectation value (n\ap\n) vanishes. Therefore, the 
one-gluon exchange kernel can only arises from the second term in Eq. (226) where the Green function Q{XTp; S^rja; E) 
reduces to the ordinary one G{Xtp; ^rja; E) in the lowest order approximation because the last term in Eq. (226) 
vanishes in this case for the reason as argued for the function i?(Arp; 7+(5~). To evaluate the inverse S~^{iiv; 7+5~), 
we first evaluate S{'-^^5~ ; iiv). Prom the expression denoted in Eq. (179), it is easy to find that the nonvanishing 
contribution of S{'y^5~;fjbu) is given by 



(243) 



where S^+i,+ and ^5-^- are the quark and antiquark equal- time propagators. Their expressions, according to the 
common definition, can be read from Eqs. (139) and (140) by setting n — T2 — > 0+, 



Therefore, we can write 



where 



Ss-iJ.- =Ss-ij,-\{q2,0^) = Ss-^^-l[nf{q2)-nf{q2)]. 



S-\iJi-u+;^+5-) = 6.^S^-s-S-\j6) (245) 



5(7^) = 1 - i [n/ (^1 ) - n/ (gi ) + n/ (0*2 ) - n/ (5*2 )] (246) 

Thus, to derive the lowest order approximate kernel, we only need to consider 

K{a+(3-;j+d-;E) = ^ ^ /(a+r ; ATp)G(Arp; /x^^a; i;)/(Mi^a; (5-7+)5-i(7<5). (247) 
Xrp {Co- 
Prom Eqs. (169), (207), (80) and (67), it is clearly seen that when A, r, /z and u take A+, t~ , /i~ and the functions 
f{a^P~', Xrp) and f{pva; 5~ ^'^) give the quark-antiquark interaction taking place in the t-channcl scattering process; 
while, when A, r, /x and u take A~, r+, fi^ and v~ , the /(a+/3~; Arp) and f{p,va] 5~7^) will give the quark-antiquark 
vertices which describe the qq annihilation process. Since the expectation value of the qq color matrix appearing 
in the qq lowest order annihilation process is zero in the color singlet, it is only necessary to consider the following 
interaction kernel 

X(a+r;7+<5-;£^) = E E Y.[f{^^p-:>^^r-p+)G{X+T-p+-iJi-v+a--E)f{p-v+a--6-^+) 

tiu\Tpir,a (248) 

+f{a+(3-;X+T- p-)G{X+T- p-;ii-u+a+;E)f{p.-u+a+;5-^+)]S-\^5). 

Noting that the fimctions f{a^f3~; X'^t~ p) and f{p~v^a\ 5~^'^) are proportional to the coupling constant g. in the 
approximation of order , we only need to consider the zero-order of the Green function G{X'^t~ p^\ijrv^a^]E). 
This Green function may easily be derived from the generating functional Z^\j] represented in Eqs. (96), (119), (138) 
and (150). The result is 

^/,+ __ - + + _ \_ J S'^Z°\j} I 

= 6x^Aq{kx,n - T2)Sr^Aq{kr,Ti - T2)d Ag{kp, Ti - T2) (^^9) 

= G{X+T~ p+\p,~v+a~\Ti - T2) 

The last equality in the above arises from the fact that Ag{kp,Ti — T2) is an even function of ti — T2. With the 
expressions given in Eqs. (124) and (140), the Fourier transform of G(A+r~p~; //"i^+a^; n — T2) can be found to be 
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G(A+T p ;m u+a+;uJn)=5\u5r^Sp„^/:i.{XTp) (250) 

where A(Atp) will be specified soon later. 

Substitution of Eqs. (169), (207) and (250) in Eq. (248) leads to 

K{a+(3-;^+5-;E) = A{a+X+ p+)A{X+^+ p-)I^{\(3p)5ps 

+ Y.MT-p-p+)A{5-T-p~)I\{aTp)5,,, ~^A{a+-i+p+)A{5-p-p-)^{^pp) (251) 

-Y.A{5-l3-p+)A{a+-i+p-)l^{a6p)]S^^{-i5). 
p 

where the first two terms on the right hand side are unconnected and represent the quark and antiquark self-energies, 
while the remaining two terms precisely give the t-channel one-gluon exchange kernel. In view of the expression in 
Eq. (80) and the definitions in Eqs. (67) and (71), we can write 

A(a+7+p±) = ig{2T^f5\p^ - q, T k)u., {pi)T'w^ {q^){2T,)-^'\2u{k))-^'^el{k) (252) 

and 

A{5-p-p^) = ig{2^f3H<h -P2 ± fc)F,^(g2)r-7M^'..(P2)(27r)-3/2(2c.(fc))-V2eA(^) 

= -igi27rfS^{q2^P2 ± fc)«.,(|g)(-T«)7.«.^(92)(27r)-3/2(2^(fc))-V2e^'(fc) ^ ' 

where the last equality in Eq. (253) is obtained by the charge conjugation transformation. In the above, we have set 
ka = Pi) = P2, — qi, ks = q2 and kp = k. Prom Eqs. (249), (250), (124) and (140), it is easy to get 

A{aSp) = i/f,drp'-"-A„(pi,r)A,(Q2.r)A,(A^.r) 



illill,—. — ILrltcIbi — ■ \ 

J J O tUJ^—Sa — SS—'-^p J J O lUJn—Sa—SS-\-^p 

J J li^„ — Ea+£S—'^p J J tlAJn—Ca+es+iAJp 

—nhifn',-. ] + nfTifTV^-. ; ^ — - 



(254) 



and 



A(pi,g2, k;ujn) 



A(7/3p) = i/f^dTe*-"-A,(p2,T)A,(9-i,r)A,(fc,T) 

J i b iu>n-e~,-ei3-Up J J b lujj^-e^-ep+tjj p 

= A{p2,qi,k;Wn) 



where 



= y^pf +mi, 50 = y/pf +rn2, = y^c/f +mi, eg = +m2, Up = k (256) 



It is noted here that the chemical potential is not taken into account for the bound state. Another expressions of the 
functions A{a6p) and A{'yPp) may be given by making use of the expansions presented in Eqs. (AlO) and (A19) in 
the Appendix. Since the expressions contain infinite series, it might be not convenient for our purpose. Upon inserting 
Eqs. (252)- (255) into the last two terms in Eq. (251) and noticing the definition in Eq. (72), after completing the 
integration over k and the summation over the polarization index, we obtain 

K{a+l3-;j+S-;E) = K{pi,p2]qi,q2]E) 
= {2n)-^5'^{pi+P2- qi- q2)V{pi,p2;q\,q2\E) ^ ' 

where the ia;„ in Eqs. (254) and (255) has been replaced by E and 
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V{pi,P2;qi,q2;E) = Ua^{pi)T''^^u„'^{qi)ua^{p2)T''*^^u„i^{q2)D{pi^p2\ (258) 

in which 

D{pi,P2;qi,q2) = „ — ir^[A(pi,g2,Pi - gi) + A(p2,gi,Pi - qi)]S^^{qi,q2) (259) 
ou)[pi — qi) 

With the kernel given above, the equation in Eq. (192) can be written as 

[E-e{pi) - e{p2)]xPaiPi,P2) = j d^qid^q2K{pi,p2;qi,q2)xPaiqi,q2) (260) 



When we introduce the cluster momenta 

P = Pi+P2, Q = qi+q2, q ^ V2P1 - 'niP2, 
k = 77291 - r?ig2, VI = V2 = 

in the center of mass frame, Eq. (260) will be represented as 



(261) 



^y(P,g,fc)xPa(fc). (262) 

As we know, the spinor function in Eq. (258) can be written as u^(p) =^ U{p)(pa where (p^ is the spin wave function 
and U {p) is the ordinary Dirac spinor determined by the Dirac equation. The spin wave functions may be absorbed 
into the amplitude ■^(q:/3; E) appearing in Eq. (193). Thus, corresponding to the equation in Eq. (262), the equation 
in Eq. (193) will be represented as 

[E - ei(q) - £2(g)]V'Pa(«) = j ^(^; k)i^Pa{k) (263) 
where ipPa{q) stands for the color singlet wave function given in the Pauli spinor space and 

V{P; g, k) = -^9^U{pi)-i^,U{qi)U{p2)rU{q2)D{P, q, k) (264) 

represents the one-gluon exchange interaction Hamiltonian which formally is the same as given in the case of zero- 
temperature [17]. In Eq. (264), we have recovered the Minkowski metric for the 7— matrix in order to compare with 
the ordinary zero-temperature result and considered that the expectation value of the gg'color operator Tf (— T2**) in 
the color singlet equals to — |. 



VIII. CONCLUDING REMARKS 



In this paper, there arc two new achievements. One is that the path-integral formalism of the thermal QCD has 
been correctly established in the coherent-state representation. The expression of the QCD generating functional 
formulated in this paper not only gives an alternative quantization of the QCD, but also provides a general method of 
calculating the partition function, the thermal Green functions and thereby other statistical quantities of QCD in the 
coherent-state representation. In particular, the generating functional enables us to carry out analytical calculations 
without being concerned with its discretized form. As one has seen from Sect. 4, the analytical calculation of the zero- 
order generating functional is more simple and direct than the previous calculations performed in the discretized form 
given either in the coherent-state representation or in the position space [18-22]. The coherent-state path-integral 
formalism corresponds to the operator formalism represented in terms of creation and annihilation operators. In 
comparison with the latter formalism which was frequently applied in the many-body theory [21, 29], the coherent- 
state path-integral formalism has a prominent advantage that in calculations within this formalism, use of the operator 
commutators and the Wick theorem is completely avoided. Therefore, it is more convenient for practical applications. 
It should be noted that although the QCD generating functional is derived in the Feynman gauge, the result is 
exact. This is because QCD is a gauge-independent theory. As shown in Sect. 4, in the partition function derived 
in the Feynman gauge, the unphysical part of the partition function given by the unphysical degrees of freedom of 
gluons is completely cancelled out by the partition function arising from the ghost particles. Certainly, the generating 
functional formulated in the coherent-state representation can be established in arbitrary gauges. But, in this case. 
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the gluon propagator would have a rather complicated form due to that the longitudinal part of the propagator will 
involve the polarization vector. Another point we would like to mention is that to formulate the quantization of the 
thermal QCD in the coherent-state representation, we limit oursclf to work in the imaginary-time formalism. It is no 
doubt that the theory can also be described in the real-time formalism. We leave the discussion on this subject in the 
future. 

The main achievement of this paper is the foundation of a rigorous three-dimensional equation for the qq bound 
states at finite temperature. Especially, the interaction kernel in the equation is given a closed expression in the 
coherent-state representation. This kernel, as shown in Eq. (226), contains only a few types of Green functions with 
some definite coefficients. We also give the corresponding closed expression represented in the position space. As shown 
in Eqs. (231)-(242), in this expression only a few types of Green functions and commutators are involved. Therefore, 
the kernel can not only easily be calculated by the perturbation method, but also is suitable for nonperturbative 
investigations by using the lattice gauge approach and some others. Since the kernel contains all the interactions 
taking place in the bound state, obviously, the kernel and the equation presented in the preceding sections are much 
suitable to study the quark deconfinement at high temperature which is nowadays an important theoretical problem 
in the high energy physics. It is expected that an accurate nonperturbative calculation of the kernel could come up 
in the future so as to give the problem of quark deconfinement a definitet solution 
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X. APPENDIX: DERIVATION OF THE GENERATING FUNCTIONAL REPRESENTED IN THE 

POSITION SPACE 

To confirm the correctness of the results derived in Sect. 4, in this appendix, we plan to derive the familiar pertur- 
bative expansion of the thermal QCD generating functional represented in the position space from the corresponding 
one given in section 4. For this purpose, we need to derive the generating functional represented in position space for 
the free system which can be written as 

z'[J\ = Zl[j;]Zl[lJ]Z^^[K\K''] (Al) 

where Zg[J^], Z^[I,I] and Z^[K"-,TC°'] are the position space generating functionals arising respectively from the 
free gluons, quarks and ghost particles and J^, I, I, if" and if" are the sources coupled to gluon, quark and ghost 

particle fields respectively. In order to write out the Z^[I, I], Zg[J^] and Z"[if if"] from the generating functionals 
given in Eqs. (119), (138) and (150), it is necessary to establish relations between the sources introduced in the 
position space and in the coherent-state representation. Let us separately discuss the functionals Zg[I,I], Zg[J^] and 
Z°[K'',K ]. First we focus our attention on the functioinal Z^ [I, I] . Usually, the external source terms of fermions 

in the generating functional given in the position space are of the form Jq dr J d^x[I{x,T)ip{x,T) + 'iI){x,t)I{x,t)] 
[22, 32]. Substituting in this expression the Fourier expansions in Eqs. (58) and (59) for the quark fields and in the 
following for the sources 



^i(f,r)=/(5f^^/(fc,r)e^'=- 

Hx,T)=Jj0j^I{k,T)e 



("I^ - ' -E- (A2) 
we have 

/(f dr J d^x[I{x, r)V'(x, r) + V^(x,j)i(f , r)] drj d^k[r]* (fc, T)bs{k, r) ^^^^ 



+bl{k,T)r]s{k,T) + T]l{k,T)ds{k,T) + dl{k,T)r]^{k,T)] 
where 

■qs{,T) =Jis{k)I{k,T), ■r]*{k,T) = 7{k,T)us{k), 
r]Ak,T) = -I{-k,T)Vs{k),rjt{k,T) = -v,ik)Ii-k,T). ^ ' 

From these relations and the property of Dirac spinors, the relation in Eq. (137) is easily proved. On substituting 
Eq. (A4) into Eq. (138), one can get 
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rP r0 r _ 

Z^[I, I] = exp{ / fin / dT2 / d^kl{k, n)SFik, Ti - r2)/(fc, T2)} (A5) 
Jo Jo J 

where 

Spik, n - T2) = [(k + m)/e{k)]Ag{k, n - T2) (A6) 

with k =7(ufc'' and k^ = {k, iSn)- By making use of the inverse transformation of Eq. (A2), the generating functional 
in Eq. (A5) is finally represented as [22, 32] 

Z°[I,I]= Z°exp{ ct^xi d'^X2l{xi)SF{xi-X2)I{x2)} (A7) 
Jo Jo 

where x = {x,t), d'^x = drd^x and 

d^k 



Sf{xi -X2)= j (^^^(^'^1 - (A8) 

It is well-known that the propagator A** (fc, n — T2) is antiperiodic, 

{k, n - T2) = - A-' (fc, n - T2 - if n >■ T2 

Af {k, n - T2) = -Af (k, Ti-T2+(3), if Ti^T2 ^ ' 

This can easily be proved from its representation in the operator formalism as shown in Eq. (46) or Eq. (154) with 
the help of the translation transformation 6s(t) = e'^^ bse~^^ and the trace property Tr{AB) = Tr{BA). According 
to the antiperiodic property of the propagator, we have the following expansion 

^g{k,T) = ^^A,{k,en)e-''''^ (AlO) 

n 

where r = ri — T2, e„ = ^(2n + 1) with n being the integer and 



r/3 



e{k) 



A,{k,en) = I dre^'-^ A,{k,T) = —^^^ (All) 

Jo £j^ + e[kY 

Substituting Eq. (AlO) into Eq. (A6) and noticing the above expression, it can be found that [22, 32] 

I r (^k e'^-(^i-^2)-fe„(ri-r2) 

Sf{xi -x2) = -ry / 7—^ = A12 

which is the familiar expression of the thermal fermion propagator in the position space. 

Next, we discuss the generating functional [J°]. The source term of gluons in the generating functional given in 

the position space is commonly taken as dr J d^xJ'^{x, t)A"''^{x, t) [22, 32]. Employing the expansions in Eq. (60) 
and in the following 

/d^k 
^^-p7^j;(fc,T)e*-, (A13) 

we can write 

J^' dr j d'xr^ix, t)A'^{x, t) = l^drj d'k[a*{k, T)al{k, r) + a'^^{k, r^k, r)] (A14) 

where 

C,{k,r) = (2a;(fc))-i/2e^(^)j;(fc,r) = in-k,r) (A15) 
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in which the last cquahty foUows from that the J^{x,t) is a real function. Inserting the relations in Eq. (A15) and 
then the inverse transformation of Eq. (A13) into Eq. (119) and considering completeness of the polarization vectors, 
one may find the generating functional Zg[J^] such that [22, 32] 



Z°[J^] = Z° exp{i J^^ d^x^ j^^ d^x^j;{x^)Df,{xr - X2)Jt{x2)} (A16) 



d^x^ / 
/o Jo 

where 

/^Sju 1 _^ 
^-^A,(fc,n - r2)e*(---) (A17) 

By the same argument as mentioned for Eq. (A9), it can be proved that the gluon propagator Ag{k,Ti — T2) is a 
periodic function 



Ag(fc, n - T2) = Ag(fc, T1-T2- /3), if Ti >- T2; 
Ag{k, n - T2) = Ag{k, n - T2 + (3), if n -< t2. 



(A18) 



Therefore, we have the expansion 

^9ik,r) = i^A,(fc,u;„)e---- (A19) 



where a;„ = and 



Ag{k,u;n)= [ dre-"-A,(ft,r) = -^^^^. (A20) 
Jo s^+Lo{ky 

Upon substituting Eqs. (A19) and (A20) in Eq. (A17), we arrive at [22,32] 

D,AX, -2)-% (2^)3^2+^(^)2^ (^21) 

which just is the gluon propagator given in the position sj>ace and in the Feynman gauge. 

Finally, wc turn to the generating functional [K'^ ,K ]. In accordance with the expansions in Eqs. (62) and (63) 
for the ghost particle fields and those for the external sources: 

i^"(5,r) = /(5^K"(fc,r)e*-^ 

K\S,r) = ! j0j,K\k,r)e-^^-^ ^""''^ 

the relation between the sources in the position space and in the coherent-state representation can be found to be 

/(f rfT/d3a;[i?"(f,r)Cnf,T) + C"(a-,T)/^(f,r)] = dr iyk[C,:{k,T)ca{k,T) 
+cl{k, r)Uk, t) + Cik, T)ca{k, t) + cl{k, r)C,{k, r)] 

where 

Uk,r) = (2u;(fc))-V_2i^-(A,r), Q{k,r) = (2u;(fc))- V27?"(A, r), 
Uk, t) = -{2w{k))-'I^K\-k, r), C:(^, r) = -(2a;(fc))- V2/^«(_fc, r). ^ ' 

from which the relations denoted in Eq. (149) directly follows. When the above relations and the inverse transforma- 
tions of Eq. (A22) are inserted into Eq. (150), one can get 

Z°[K'',k"] = ZOexp{- Td^xi f d'x2K''{xi)Af{x^ - X2)K\x2)} (A25) 
^0 Jo 

where 
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l2¥F^,^^^^-j,- 



Ia6J_Y^ f d^fc 1 pik-(xi—X2)—iu)n(T-L — T2) 



(A26) 



which just is the free ghost particle propagator given in the position space [22, 32]. fn the last equality of Eq. (A26), 
the expansion given in Eqs. (A19) and (A20) have been used. 

With the generating functionals given in Eqs. (A7), (A16) and (A25), the zeroth-order generating functional in 
Eq. (Al) is explicitly represented in terms of the propagators and external sources. Clearly, the exact generating 
functional can immediately be written out from Eq. (95) as shown in the following 

Z[J] = exp{- d^xnj{j^^)}ZV] (A27) 

where J stands for J, 7, J« K'^ and iT, represents the differentials -^M' i^fw ' 5^%) ^'^'^ ~Ik^ 

and T^i{jj(^) can be written out from Eq. (57) when the field functions in Eq. (57) are replaced by the differentials 
with respect to the corresponding sources. 
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